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Adiabatic quantum computing has recently been used to factor 56153 [Dattani & Bryans, 
arXiv:1411.6758| at room temperature, which is orders of magnitude larger than any number at¬ 
tempted yet using Shor's algorithm (circuit-based quantum computation). However, this number 
is still vastly smaller than RSA-768 which is the largest number factored thus far on a classical com¬ 
puter. We address a major issue arising in the scaling of adiabatic quantum factorization to much 
larger numbers. Namely, the existence of many 4-qubit, 3-qubit and 2-qubit interactions in the Hamil¬ 
tonians. We showcase our method on various examples, one of which shows that we can remove 
94% of the 4-qubit interactions and 83% of the 3-qubit interactions in the factorization of a 25-digit 
number with almost no effort, without adding any auxiliary qubits. Our method is not limited to 
quantum factoring. Its importance extends to the wider field of discrete optimization. Any CSP 
(constraint-satisfiability problem), psuedo-boolean optimization problem, or QUBO (quadratic un¬ 
constrained Boolean optimization) problem can in principle benefit from the "deduction-reduction" 
method which we introduce in this paper. We provide an open source code which takes in a Hamilto¬ 
nian (or a discrete discrete fimction which needs to be optimized), and returns a Hamiltonian that has 
the same unique ground state(s), no new auxiliary variables, and as few multi-qubit (multi-variable) 
terms as possible with deduc-reduc. 

PACS numbers: 06.20.Jr, 31.30.jh, 31.50.Bc , 95.30.Ky 


I. INTRODUCTION 

The quantum algorithm which has perhaps gener¬ 
ated the most enthusiasm about quantum computing, 
is Shor's algorithm for factoring integers Eli. How¬ 
ever, despite being celebrated for more than 20 years, 
this algorithm has still never been successfully imple¬ 
mented for determining the factors of any integers with¬ 
out using knowledge of the answer to the problem Hi. 
By using knowledge of the answer to the factorization 
problem, one can choose a base such that Shor's algo¬ 
rithm can be implemented with fewer qubits, and by 
doing this, the algorithm has successfully been imple¬ 
mented for factoring 15 M and 21 Iia. However, at 
least 8 qubits are needed for genuinely factoring 15 with 
Shor's algorithm a, and the largest number of qubits 
ever successfully used in the algorithm was 7 H- 

Adiabatic quantum computing (AQC) has succeeded 
in factoring much larger numbers with far fewer qubits, 
without any assumptions about the answer to the fac¬ 
torization problem 0]. The largest number found so far 
that has been factored by the room-temperature AQC 
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discrete minimization e)weriment of lITlIl is 56153 and 
it only needed 4 qubits IJ]. Furthermore, quantum an¬ 
nealing can be used on AQC algorithms with up to 2048 
qubits, which is enough to factor the 100-digit number 
RSA-100. 

Furthermore, quantum armealers have doubled in 
number of qubits every year between the 4-qubit Ca¬ 
lypso machine in 2005 and the 2048-qubit Washington 
machine in 2015 (a phenomenon analogous to Moore's 
law and known as "Rose's Law"). RSA-220, which is 
the smallest RSA number that has not yet been factored 
by any computer (whether quantum or classical), would 
need only ~ 5000 qubits to factor successfully with the 
AQC algorithm of |0], which would be well within the 
8192-qubit capacity of a 2017 quantum annealer if Rose's 
law continues to hold as it has done for the last 10 years. 

However, there is a major obstacle holding back such 
quantum annealers with colossal numbers of qubits, 
from implementing the AQC factoring algorithm. All 
such quantum annealing devices reported to date can 
only implement AQC algorithms which have at most 
2-qubit interactions in their Hamiltonian. The algo¬ 
rithm which can factor RSA-220 with ~ 5000 qubits 
has a Hamiltonian with many 4-qubit and 3-qubit in¬ 
teractions. In 2007 Schaller and Schutzhold devised an 
alternate AQC algorithm for factoring integers, which 
only has up to 2-qubit interactions in the Hamiltonian, 
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and whose Hamiltonian's spectral width (which is used 
for determining the runtime of the annealing) increases 
only pol 5 momially in the number being factored llT3l . 

The problem with this AQC algorithm that uses at 
most 2-qubit interactions and whose spectral width only 
increases polynomially, is that it requires adding a large 
number of auxiliary qubits (27225 for RSA-100, and 
133 225 for RSA-220). In this paper we present a method, 
called "deduc-reduc" for reducing fhe number of 4- 
qubif and 3-qubit terms in the AQC Hamiltonian without 
adding auxiliary qubits. 

While the motivation for this work was the problem 
of factoring numbers using AQC, deduc-reduc can be 
used for a wider range of discrete optimization prob¬ 
lems such as constraint satisfaction problems (CSPs) and 
pseudo-Boolean optimization problems, and problems 
arising in the field of artificial infelligence/neural nef- 
works Cl. All of these problems can benefit from quan¬ 
tum armealing because of the capability for quadratic 
speed-up over simulated annealing m. Quantum an¬ 
nealing has also recently been shown to be able to pro¬ 
vide exponential speed-up for some specific discrete op¬ 
timization problems Cl/ the potential for speeding 
up more general discrete optimization problems is an 
active area of research. 

It has also been shown that AQC can simulate any 
circuit-based qua ntum algorithm with only polynomial 
overhead 11^,11^ . Therefore, apart from fhe exponen¬ 
tial speed-up offered for many classes of problems, AQC 
can provide polynomial speed-up over the best classical 
3-SAT algorithms [? ] and other CSPs Cl - arid can pro¬ 
vide exponential speed-up for combinatorial optimiza- 
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tion problems for which an approximate solution is suf¬ 
ficient m . 


II. A QUICK EXAMPLE 

It has been shown with concrete examples in many 
papers (UnLCIlIi that integer factorization is equiva¬ 
lent to solving a system of binary equations such as: 


*1 + *2 + *3 = 1 

( 1 ) 

*1*4 + * 2*5 = *3 

( 2 ) 

*1 + 2*2 = *3 + 2*4 

(3) 


which is equivalent to finding the minimum of an objec¬ 
tive function formed by subtracting the right hand side 
of each equation from its left, squaring, and summing 
(and remembering that the square of a binary variable is 
itself): 

H = (a:i+a:2+a:3-l)^ + (a:ia:4+a:22:5-a:3)^+ (4) 

( 0 : 1 + 22 : 2 - 2 : 3 - 2 * 4 )^ (5) 

= 2*12:22:42:5-22:12:32:4-22:22:32:5-22:22:3+62:12:2 (6) 

- 3 * 1 * 4 - 8 * 2 * 4 +* 2 * 5 + 3 * 2 + 4 * 3 * 4 +* 3 + 4*4 + 1 (7) 

If we represent each variable Xi as a qubit based on a 
Pauli spin matrix in its own Hilbert space: 



and each term of Eq. [6]using left Kronecker products (or 
"tensor products") in the Hilbert space of all qubits, for 
example: 


o o o/O 0\ (0 0\ /I 0\ /O 0\ (0 0\ 

2*1*22:42:5 = 2*1*2132:42:5 = 21g i|®loi|®loij®loij®loir 

= diag(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,0,0,0,2), (10) 

we see quickly fhat the objective frmction of Eq. [6]can be represented as a diagonal 32x32 Hamiltonian matrix: 


Tf = diag(l, 1,5,5,2,2,10,10,4,5,0,1,3,2,3,2,1,1,2,2,2,2,5,5,10,11,3,6,9,8,4,5). 

The lowest eigenvalue is clearly 0, and it is at position: 

diag(0,0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0) = |*i = 0,*2 = 1,*3 = 0,*4 = 1,*5 = 0). (11) 


This may seem like an rmnecessary extra complica¬ 
tion to find the minimum, especially since if we have 
N variables, it involves representing all 2^ values of H 
and then searching for the lowest one. However, since 
the lowest eigenvalue of a Hamiltonian corresponds to 
its ground state, the minimum of a 2^ variable system 
can be found simply by connecting N qubits together 


with appropriate energies and coupling strengths, then 
finding fhe ground state of fhis system. This avoids hav¬ 
ing to evaluate the function 2" times. But it is much 
harder to couple three qubits together as in the cubic 
term 2 *i* 2 * 3 , than it is to couple two qubits together 
as in the quadratic terms, which is of course sfill harder 
fhan not having to enforce any couplings (buf if only 
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linear terms exist, the problem can be solved trivially 
by setting all negative terms to 1 and positive terms to 
0). Likewise, even if solving this problem on a classical 
computer, multi-variable (or multi-qubit) terms make 
the problem much more difficult, and it is desirable to 
eliminate any terms containing more than 2 variables, 
since this would transform the problem into a quadratic 
unconstrained Boolean optimization (QUBO) problem, 
whereby we can use many beautiful algorithms and re¬ 
sults such as the fact that the solution can be found in 
polynomial fime if all quadratic terms are negative 10 - 
124)1 ■ Furthermore, it is also desirable to eliminate as 
many quadratic terms as possible. 

A. The method 

We now demonstrate how deduc-reduc can be used to 
eliminate multi-qubit interactions (or high-order terms). 

From Eq. [T]alone we can make the deductions: 

XiX2=X2X3 = X3Xi=0 ( 12 ) 

since the LHS of Eq. [T]cannof be greater than 1, and this 
would not be true if any of fhe products in Eq. [121 were 
1 . 


Naive substitution 

Suppose that we substituted Eq. [T2]back into H: 

H = -‘iXlXi-9>X2Xi+X2X5 + ‘iX2+'^XzXi + X3+AXA + l. (13) 

Then we find thaf while iT(0,1,0,1,0) = 0 matches the 
original Hamiltonian as required, we have also intro¬ 
duced a new state with H = -3 and another for which 
H = -2. This is a problem, because if the minimum is 
not 0 , then the ground state does not encode the solution 
to Eqs[T][3l which is the integer factorization problem we 
want to solve. 


Deduc-reduc 

Instead, we can look at each term we would like to 
substitute Eq. [12] into, and calculate an associated error 
term. 

We want an error term such that it is equal to 0 if 
xiX 2 = 0 and strictly greater than 0 ifxia ::2 + 0. In this 
particular example, we can take xiX 2 to be our error 
term and similarly for X 2 X^ and X 3 ,xi. We will now work 
through the rest of this example in order to illustrate the 
method, but a general formulation can be formd in Sec- 
tionUni 

We have three cubic/quartic terms to consider: 

1. 2 xiX 2 XiX^. Here we make use of fhe judgment 
xiX 2 = 0. Eor any state, 2 xiX 2 XiXr, can only evalu¬ 
ate to 0 or 2. In the former case, a straight substitu¬ 
tion does not make any difference to the value of 


H and sefting fhe term to 0 is permissible. If the 
term evaluates to 2 (i.e. all variables are equal to 
1 ) and we perform a sfraighf substifution, we are 
taking 2 away from the energy of the state. This 
may introduce new states with iT = 0 or even neg¬ 
ative valued states. To compensate for fhis, we add 
2 a;ia ;2 back on fo fhe original Hamiltonian to pre¬ 
serve the original value of H. We have then made 
the transformation 2 X 1 X 2 X 4 X 5 -»■ 2 a;ia; 2 . While this 
step may add 2 to the original energy of a sfate, 
the important detail is that we know this will not 
happen for any ground stafes. 

2. -2X1X3X4. This can only fake the values 0 or -2. In 
the first case, a straight substitution makes no dif¬ 
ference. In fhe second, a sfraighf substifution for 
0 would add 2 to the value of the incorrect state, 
which is safe to do, (it is even desirable, since in¬ 
creasing the energy gap between the ground and 
excited states can allow the ground state to be 
found quicker by armealing or adiabafic evolufion 
methods, and the power to manipulate the energy 
landscape like this is the subject of our forthcom¬ 
ing paper lH^). Hence we can make the transfor¬ 
mation - 2 x 12 : 33:4 -» 0 . 

3. - 2 X 2 X 3 X 5 . As with case 2, we make the transfor¬ 
mation - 2 X 2 X 3 X 5 -> 0 . 

Bringing the above substitutions together we get a 
Hamiltonian with the same number of qubits as before, 
but has no 4-qubit (quartic) nor 3-qubit (cubic) terms, 
and whose only ground state is precisely the ground 
state of the original Hamiltonian: 

H = 8xiX2-2X2X3-3xiX4-SX2X4+X2X5+ (14) 

3x2+ 4 x 32 : 4 + 2 : 3 +4x4+1. (15) 

A quick numerical check shows that our new Hamil¬ 
tonian has a single zero at exactly the original solution 
as required. We have therefore used our deductions in 
Eq. |T2|to do a reduction of fhe 4-qubif and 3-qubif terms. 


Removing quadratic terms 

Notice that when the coefficient of a term was nega¬ 
tive, we were able to perform a straight substitution, as 
doing so only added to the energy of a state and no er¬ 
ror term was needed. This can also be used to eliminate 
quadratic terms like - 2 x 2 x 3 , to provide extra simplifica¬ 
tion. Our final Hamiltonian is then: 

ff = 8x1X2 -3xiX4 -8x2X4+X2X5+3x2 +4X3X4 + X3+4X4 + 1. (16) 

Note that in general we can remove quadratic terms by fol¬ 
lowing the method described in Section |nl| However, this is 
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Table I. The three examples used to illustrate the power of 
deduc-reduc for removing multi-qubit interactions. 


Example 

Product 

Length of product 
in binary 

1 

455937533473 

40 

2 

292951160076082381 

60 

3 

1208925727750433490141601 

80 


difficult in practice because any deduction must be at least 
quadratic (else we perform a straight substitution) and re¬ 
quires an appropriate error term. 


B. Performance of deduc-reduc on large examples 

We now consider the performance of deduc-reduc on 
Hamiltonians whose ground states encode the prime 
factors to larger numbers, displayed in Table U 

For each example, we generate the carry equations 
that express the factorization and then form fhe Hamil¬ 
tonian as in SHi: this is our starting Hamiltonian Hq. 
Next, we perform simple logical deductions to solve for 
the most obvious variables (as in Eq. [121 and the Supple¬ 
mentary Material of lITlIl ). using no enumeration of so¬ 
lutions or deduc-reduc machinery. We call this reduced 
Hamiltonian Hi. 

Finally we perform searches through the state space 
by enumerating plausible solutions and looking for pat¬ 
terns that can be used by deduc-reduc. We perform a 
breadth-first search that halts when it has found more 
than n plausible solutions. We call the Hamiltonian 
formed using deduc-reduc with the deductions from 
such a search iJ„. 

Enumerating all possible solutions is akin to doing 
a brute force search for the minimum, so of course it 
would make no sense to run this algorithm for n any¬ 
where close to the total size of the search space 2 ^ where 
N is the number of qubits in the Hamiltonian. How¬ 
ever, in the examples from Table HI we show in Table HI] 
that with n vastly smaller than 2^, we can still eliminate 
a significant number of multi-qubit interactions Q. No 
deduc-reduc example we present takes more than a few 
hours of CPU time on an average laptop, and 1000 states 
takes just a few seconds. The code is written in Python 
with scope for optimization and parallelization. 


^ While the factorization problem can be directly solved by brute force 
using 2^" < 2^ trials, where m is the bit-length of the number we 
wish to factorize, most discrete optimization problems will not have 
this feature, so the brute force search to solve the problem will gen¬ 
erally require 2^ trials where N is the number of variables. Further¬ 
more, the number of states in our search is still vastly fewer than 
even 2™-. 


In Example 1 we are unable to make a big difference 
in the number of multi-qubit interactions because both 
the simple judgments and local searches are imable to 
find many solutions or deductions that are useful for 
deduc-reduc. In Example 2, the large improvement in 
the qubit profile comes from the power of simple judg¬ 
ments, which leaves little 'easy' work for the search and 
deduc-reduc to do afterwards. The remarkable success 
of deduc-reduc in Example 3 stems from the effective¬ 
ness of the breadth-first search's ability to find patterns. 
The factors that affect the performance of the simple 
judgments and the breadth-first search are a potential 
area for future research though we believe it is related 
to the distribution of Os and Is in the binary representa¬ 
tions of the factors. 


III. GENERAL TREATMENT 

A. General formulation 

Let ® be a binary string of length n. Suppose we have 
a general Hamiltonian H{x) and a deduction f{x) = 
g{x) that is true for all ground states, and f,g are poly¬ 
nomials with the degree of g strictly less than the de¬ 
gree of /. Suppose futher that we are able to write 
= g{x)f{!>^ + r{x) for some polynomials q + 0 and 
r. Then substitution of / for g would give a new poly¬ 
nomial with fewer high order terms. 

To allow us to perform this substitution and preserve 
the global minima, we add an error term XC{x) of suf¬ 
ficient magnitude to compensate for the substitution. If 
we start with the deduction "/(x) = g{x) for all minima 
x" and form a corresponding error term C = (f-g)^, we 
define the new pol 5 momiaI: 

H'{x) = q{x)g{x) + r{x) + XC{x) (17) 

= Qix)g{x)+r(x) + X(f-gf (18) 

where |( 7 (s)| < A for all states x. We claim that H' has 
exactly the same minima as H, and no new minima, and 
we prove it in the next sub-section. 

B. Proof of formulation 

To show that this H'{x) has the required properties, 
we condition on whether or not a state x is a ground 
state of the original Hamiltonian. 

1. If X is a ground state then f{x)=g{x) and C'(x) = 
{f-g )^=0 so we have 

H'{x) = q{x)g{^+r{x) + XC{x) (19) 

= Qix)f{x)+r(x) (20) 

= H{2L) ( 21 ) 

as required. 
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Table II. Performance of deduc-reduc on the three examples from Table|I] 


455937533473 = 524309 x 869597 = 10000000000000010101 x 11010100010011011101 


Hamiltonian 

# of qubits # of deductions 

4-qubit terms 

3-qubit terms 

2-qubit terms 

1-qubit terms 

Original, Ho 

174 

1785 

3318 

1783 

150 

Simple Judgments, Hi 

148 

1750 

2915 

1407 

128 

Reduction with 100 states, TTioo 

146 88 

1645 

2828 

1404 

128 

Reduction with 1000 states, tTiooo 

144 121 

1645 

2794 

1411 

126 

Reduction with 10 000 states, Hioooo 

141 131 

1645 

2704 

1362 

123 

Reduction with 100 000 states, Hwo ooo 

138 242 

1645 

2591 

1370 

120 

292951160076082381 = 539152967 x 

543354443 =100000001000101101001001000111 x100000011000101110111001001011 

Hamiltonian 

# of qubits # of deductions 

4-qubit terms 

3-qubit terms 

2-qubit terms 

1-qubit terms 

Original, Hq 

294 

6930 

8816 

3541 

268 

Simple Judgments, Hi 

200 

3686 

5061 

2222 

180 

Reduction with 100 states, Hwo 

199 100 

3388 

4942 

2193 

180 

Reduction with 1000 states, Hwoo 

199 160 

3388 

4924 

2192 

180 

Reduction with 10 000 states, Hioooo 

189 111 

3388 

4860 

2074 

173 

Reduction with 100 000 states, Hwo ooo 

185 129 

3388 

4747 

1992 

169 


1208925727750433490141601 = 1099511555521 x 1099511616481 

1111111111111111111111101110010111000001 X1111111111111111111111111101001111100001 


Hamiltonian 

# of qubits 

# of deductions 

4-qubit terms 

3-qubit terms 

2-qubit terms 

1-qubit terms 

Original, Ho 

430 

- 

17575 

17762 

5891 

382 

Simple Judgments, Hi 

367 

- 

16133 

16005 

5256 

367 

Reduction with 100 states, Hwo 

333 

56 

8469 

14203 

6881 

327 

Reduction with 1000 states, Hwoo 

257 

289 

1449 

3649 

6180 

257 

Reduction with 10 000 states, Hwooo 

253 

285 

1449 

3565 

6106 

253 

Reduction with 100 000 states, Hwo ooo 

236 

236 

1142 

2961 

5732 

236 


2. If X is not a solution, we have that |g(x)| < A, by 
construction, and |g(x) - f{x)\ < (gix) - f{x)) , 
as we are dealing with polynomials taking integer 
values. Hence 

0 < qix)igix)-f{x)) + \{g{x)-f{x)f. (22) 

Then, 


H'ix) = q(x)g{x) + r(x) + XC(x) (23) 

= (q(x)f{x) + r{x)) + q{x){g{x)-f{x))+ (24) 

XC{x) (25) 

=H{x)+q{x)(g(x)-f{x))+X(g(x)-f(x)f{26) 
> 0 

since H{x) > 0. 

Hence H{x) = 0 if and only if H'{x) = 0. 


C. Calculation of A 

While maximizing an arbitrary integer valued poly¬ 
nomial in n variables might seem suspiciously similar 
to the original problem, we only need an upper bound. 
Such a bound can quickly be found by the triangle in¬ 
equality. However, this might lead to much larger co¬ 
efficients than are needed. It is perfectly permissible to 
split up q{x) and apply repeatedly, as we did in our ex¬ 
ample to find Hamilfonian [141 

In the most general case, A can be picked even more 
precisely depending on the /, g, C and q, as long as the 
following inequality is satisfied: 

0 ^ g{x){g{x)- f{x)) + '^C'(x) for all states x. (27) 

IV. CONCLUSION 

For any discrete optimization problem, we have 
presented a method (called "deduc-reduc") that can 
re-formulate the objective function such that without 
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adding auxiliary variables, the new objective fimction 
has fewer high-order terms. Any input that minimizes 
the original objective fvmction also minimizes the new 
one, and no new inputs can minimize the objective func¬ 
tion. This can be used to reduce a high-order objective 
function into an easier one, such as a quadratic Boolean 
optimization problem (QUBO). On our examples for in¬ 
teger factorization, we were able to remove thousands 
of quartic and cubic terms in a few seconds on a laptop. 
We provide an open source code with sample input files 
inl^. 

The method we presented can also be used to turn 
a target Hamiltonian in adiabatic quantum computa¬ 
tion (AQC), without adding auxiliary qubits, into one 
with equivalent unique ground states, but with only 
2 -qubit terms (for implementations using, for example 
SQUID-based quantum armealing devices that carmot 
yet couple more than 2 qubits |27l] ') or low enough or¬ 
der terms for other AQC-like implementations such as 
NMR-based devices which can implement 3- and 4- 
qubit interactions with some more effort Gill. 

The method we presented here can also be used 
to manipulate the energy landscape of a Hamiltonian 
without changing the unique ground states. Particu¬ 


larly, it can be used to move low-lying excited states 
(local minima of the objective function) away from the 
groimd state (global minimum of the objective fimction) 
and to adjust the Hamiltonian's spectral width (range 
of the objective function) without adding more qubits. 
This would allow for adjusting the size of the spectral 
gap between the ground and first excited states of the 
target Hamiltonian, and therefore it can be used to de¬ 
crease the runtime of the adiabatic quantum computa¬ 
tion, and to also speed-up conventional discrete opti¬ 
mization calculations such as simulated armealing. This 
is the subject of a forthcoming publication |2^ . 

Finally, this is the first paper of a 2-part series on 
techniques for reducing multi-variable terms in discrete 
optimization problems. The second method is called 
"split-reduc" and is presented in G3l with examples fo¬ 
cused on the quadratization of Hamiltonians used for 
AQC-based determination of Ramsey numbers l^lsoll . 
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